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Abstract. We discuss a "stationary 1-flog" slicing condition for the construction of 
solutions to Einstein's constraint equations. For stationary spacetimes, these initial 
data give a stationary foliation when evolved with "moving puncture" gauge conditions 
that are often used in black hole evolutions. The resulting slicing is time-independent 
and agrees with the slicing generated by being dragged along a time-like Killing vector 
of the spacetime. When these initial data are evolved with moving puncture gauge 
conditions, numerical errors arising from coordinate evolution should be minimized. 
While these properties appear very promising, suggesting that this slicing condition 
should be an attractive alternative to, for example, maximal slicing, we demonstrate 
in this paper that solutions can be constructed only for a small class of problems. For 
binary black hole initial data, in particular, it is often assumed that there exists an 
approximate helical Killing vector that generates the binary's orbit. We show that 
1-1- log slices that are stationary with respect to such a helical Killing vector cannot be 
asymptotically flat, unless the spacetime possesses an additional axial Killing vector. 



1. Introduction 

Numerical simulations of binary black hole systems have recently experienced a dramatic 
breakthrough (see [H [2l [3] as well as numerous follow-up publications). Most of these 
simulations now adopt some variation of the BSSN formulation [H [5] together with the 
moving puncture method to handle the black hole singularities. 

The original puncture method, which factors out the singular behavior of black hole 
spacetimes analytically, was originally developed in the context of initial data P [TJ E]. 
When used for dynamical simulations (e.g. [9|,[T0l[TT]), this method did not lead to long- 
term stable evolutions. The breakthrough in the recent puncture simulations is based on 
a "moving puncture" approach, in which the singular terms are no longer factored out, 
and in which the punctures are allowed to propagate through the numerical grid. With 
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a set of empirically determined coordinate conditions, and as long as no grid point ever 
encounters a singular term directly, this prescription then leads to remarkably stable 
evolutions. 

It is quite surprising that the presence of singularities does not spoil the numerical 
calculation. This issue has been clarified by Hannam et al. p ^fTSlITi] (see also [E]), who 
analyzed the geometry of moving puncture evolutions for single Schwarzschild solutions. 
Typically, such evolutions start out with a slice of constant Schwarzschild time expressed 
in isotropic coordinates as initial data. These coordinates do not penetrate the black 
hole interior, and instead cover two copies of the black hole exterior, corresponding to 
two sheets of asymptotically flat "universes" , connected by an Einstein- Rosen bridge at 
the black hole horizon. Given their appearance in an embedding diagram, these data 
are often called "wormhole" data. The singularity at isotropic radius r = 0, where the 
conformal factor ip = 1 + M/ (2r) diverges, corresponds to the asymptotically fiat end of 
the "other" universe, and is therefore a coordinate singularity only. When evolved with 
the "1+log" slicing condition for the lapse [16] and a "F-freezing" gauge condition for 
the shift [IT], the solution passes through a short phase of "dynamical" evolution, but 
then settles down into a new equilibrium solution. An embedding diagram of this new 
equilibrium solution is shown in Fig. 2 of fH] and suggests the name "trumpet" data. 
The solution terminates at a non-zero areal radius, and therefore does not encounter the 
spacetime singularity at the center of the black hole. The singularity at the isotropic 
radius r = is therefore again a coordinate singularity only. This new equilibrium is 
a stationary slice of the "1 + log" slicing condition relative to the standard timelike 
Killing vector of the Schwarzschild solution. 

At least in the case of Schwarzschild, the brief "dynamical" phase of the evolution 
refers to a coordinate evolution only, since the spacetime itself is clearly static. Even 
though the spacetime is static, it is sliced, in this case, by a nonstationary foliation. 
The resulting nontrivial evolution of the coordinates leads to a time-dependence of the 
metric coefficients, which may well introduce or enhance numerical errors in dynamical 
simulations. This observation suggests that it would be advantageous to start the 
numerical evolution not with slices of constant Schwarzschild time as initial data, but 
instead with a slicing that corresponds to the late-time equilibrium solution. Stated 
differently, we would like to identify a "stationary" slicing of the spacetime that satisfies 
the 1+log slicing condition used in moving puncture simulations. All quantities should 
then remain constant during a time evolution (up to numerical truncation error), 
potentially reducing numerical errors and noise. 

For binary black hole evolutions, moving puncture simulations have typically used 
maximally sliced puncture initial data (see, e.g., [H [H]). In the limit of infinite 
binary separation, these data again approach slices of constant Schwarzschild time 
for two individual static black holes (see [19] for a perturbative analysis of this 
limit). Presumably, these data contain several sources of error, including some nonzero 
eccentricity (see [Sni El]) and the absence of the correct gravitational wave pattern 
originating from the prior binary inspiral. Another potential source of error arises 
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from the fact that the individual black holes are not in equilibrium when evolved with 
the l+log and F-freezing conditions, even at infinite binary separation. As described 
above for single black holes, a dynamical evolution will lead to the individual black 
holes evolving to a new quasiequilibrium coordinate state. This "coordinate-dynamical" 
process may again introduce numerical error, and may contribute to the "spurious 
gravitational radiation" that is often observed at early times in binary black hole 
simulations. 

It would again be desirable to construct initial data that describe black holes in a 
coordinate system that is already in equilibrium with respect to 1+log and F-freezing 
coordinates used in moving puncture simulations. Obviously, stationary coordinate 
systems can only exist for stationary spacetimes. While inspiralling binary systems 
are not stationary, they possess an approximate helical Killing vector, and it would be 
desirable to find a "stationary 1+log" slicing with respect to this Killing vector. As we 
will discuss below, a corresponding slicing condition can be formulated quite naturally 
in the context of the conformal thin-sandwich decomposition of Einstein's constraint 
equations. Given these considerations, this "stationary 1+log" slicing condition would 
appear as an attractive alternative to maximal slicing, which is often used in the 
construction of binary black hole initial data. However, we demonstrate in this paper 
that, in the absence of an axial Killing vector, asymptotically fiat solutions to the 
resulting equations cannot exist. In particular this rules out solutions describing binary 
black holes. 

This paper is organized as follows. In Section[2]we briefly review the basic equations, 
including the 1+log slicing condition and the conformal thin-sandwich decomposition 
of Einstein's constraint equations. We then solve the latter equations, subject to the 
1+log slicing condition, for Schwarzschild black holes in Section [31 In Section H] we then 
demonstrate why such solutions cannot exist in rotating spacetimes, in the absence of an 
axial Killing vector. We stress that this nonexistence result is valid only where the only 
Killing vector is a helical one. If the Killing vector is a combination of a rotational Killing 
vector and a time-translational one, the problem disappears. We briefly summarize 
in Section [51 In Appendix [31 we also include some results from dynamical evolution 
simulations of our Schwarzschild solutions. 

2. Basic equations 

2.1. The 3+1 decomposition 

We write the spacetime metric Qab in the 3+1 form 

Qatdx^dx^ = -a^dt^ + -fijidx' + (3'dt){dx^ + (3^dt), (1) 

where a is the lapse function, /3* the shift vector, and 7ij the spatial metric induced on a 
spatial slice S of constant coordinate time t. The normal on this slice is = a~^{l, —j3'^) 
(note that in our "Fortran" convention letters a,b,c, . . . denote spacetime indices, while 
letters i,j,k,... denote spatial indices). We also define the extrinsic curvature as 
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Kij = — (l/2)£n7jj, where denotes the Lie derivative along n*^ (we adopt the "ADM" 
convention of Arnowitt, Deser and Misner |22] rather than that of [23]). Expanding this 
Lie derivative we have 

dtjij = -2aKij + Dip J + Djpi, (2) 

where Di is the covariant derivative associated with 'jij. 

Einstein's equations can then be decomposed into two sets of equations for the 
field variables jij and Kij. One of these sets is the constraint equations, namely 
the Hamiltonian constraint and the momentum constraint, which constrain the field 
variables on each spatial slice S. The second set of equations is the evolution equations, 
one of which is given by equation ([2]), which determines the time evolution of the 
fields from one spatial slice to the next. The lapse and the shift, which determine how 
the coordinates evolve from one spatial slice to the next, appear only in the evolution 
equations, and have to be chosen independently before the evolution can proceed. 

Constructing initial data entails finding solutions to the two constraint equations. 
Before this can be done, a particular decomposition of the constraint equations has to 
be chosen (see, e.g., [241 [25l [26] ioT reviews). Particularly popular decompositions are 
the conformal transverse-traceless decomposition, which provides the framework for the 
black hole puncture initial data (see [6], [71 13 HE]), and the conformal thin-sandwich 
decomposition. As we will see below, the latter is more suitable for our purposes. 



2.2. Stationary 1+log slicing 

The 1+log slicing condition, which is often employed in dynamical moving puncture 
simulations, is usually written as 

{dt - f3'di)a = -2aK (3) 

(see [IS]). Here the mean curvature K is the trace of the extrinsic curvature, K = Y'^Kij. 
We point out that the presence of the advective shift term allows us to write the left 
hand side clS db derivative along the normal vector, 

n^Vaa = = -2K. (4) 

This means that this slicing is "covariant" in the sense that it does not depend on the 
choice of the shift - for a given slice, and a given initial lapse, the resulting slicing of 
the spacetime is unique. 

We would now like to construct a stationary slicing that satisfies the slicing 
condition ([3l). Clearly we can find a stationary slicing of a spacetime only if the spacetime 
is stationary, meaning that the latter possesses a Killing vector A stationary slicing 
of the spacetime is then dragged along by We can now construct a coordinate system 
by identifying the time coordinate vector with the Killing vector i.e. 

^» = t"^ = aKu'' + (5) 

This identification provides what we call the Killing lapse ax and the Killing shift 
Given a certain initial slice, and evolving this slice with the Killing lapse and shift will 



The stationary 1+log slicing 



5 



render all quantities independent of time. In the "1 +log" slicing, the initial value of 
the lapse is a free choice, and the foliation condition determines its time evolution. To 
construct a stationary slicing, the initial lapse must be the Killing lapse. The slicing 
lapse must remain equal to the Killing lapse when dragged along the Killing vector, i.e., 
when one chooses the slicing shift to be equal to the Killing shift. However, whether or 
not this Killing lapse simultaneously satisfies the 1+log slicing condition ([3]) depends on 
the choice of our initial slice. 

For a concrete example, consider the Schwarzschild spacetime, foliated by slices of 
constant Schwarzschild time t, in isotropic coordinates. We can then identify the Killing 
lapse as 

= l + M/(2r) 

and the Killing shift as (3^ = 0. Given that the spatial metric is time-independent 
on these slices we also have Kij = = K, and we can readily verify that this slicing 
does satisfy the 1+log slicing condition ([3]). Note, however, that this ax is negative 
for r < M/2. We can just as easily construct a counter example. In Painleve- 
Gullstrand coordinates we can identify the Killing lapse as ax = 1, the Killing shift as 
= (2M/r)^/^ and the mean curvature is K = {?,/2){2M /r^y/"^ . Evidently, slices of 
constant Painleve-Gullstrand times do not satisfy ([3]), nor do slices of constant Kerr- 
Schild or ingoing Eddington-Finkelstein time. 

Having chosen the slicing lapse as as equal to the Killing lapse, and the slicing shift 
equal to the Killing shift, we need dtax = dtas = 0. The condition ([3]) then reduces to 

2aK ' ^ ' 

We will refer to this condition as "stationary 1+log slicing". 

In the following we will employ the condition ([7j) for the construction of initial data. 
As we discussed above, the constraint equations alone do not provide conditions for the 
lapse and shift, so that it is not clear how ([7]) could be employed in any meaningful 
way. The conformal thin-sandwich decomposition, however, is constructed by imposing 
conditions on the time evolution of the spatial metric, which automatically introduces 
the lapse and shift. This suggests that we should adopt the conformal thin-sandwich 
decomposition to solve the constraint equations in the context of the stationary 1+log 
slicing condition. 

As an aside - which will provide a useful analogy later in this paper - we note that 
some authors have considered the condition ([3]) without the advective term, 

dta = -2aK (8) 

(which is no longer a "covariant" slicing condition in the sense that we discussed 
above). In this case, the Killing lapse associated with a stationary slicing satisfies 
this condition only if the slices are maximal, i.e. if the mean curvature vanishes K = 0. 
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For Schwarzschild, all maximal slicings can be parametrized by a parameter C (see 
[271 [281 [29] ) . The lapse, for example, is given by 

where R is an areal radius (in contrast to the isotropic radius r in ([6])). For C = we 
recover the shces of constant Schwarzschild time t, as in ([6|). As discussed by |30j, we 
could also parametrize this family of slicings by the value of the lapse on the horizon, 
oah- For C = 0, for example, we evidently have oah = 0. 

Dynamical simulations of Schwarzschild spacetimes, adopting the condition ([8|) and 
a F-freezing condition for the shift, settle down not to an arbitrary maximal slice, but 
instead to that with C = 3a/3M^/4, which has a limiting surface at R = 3M/2 (see 
|13j). To construct initial data that lead to a trivial time-evolution under the condition 
([8]), we would therefore use this particular maximal slice (see [31]). We will return to 
this issue in Section [31 

2.3. The conformal thin- sandwich decomposition 

We begin by conformally decomposing the spatial metric as 

%j = (10) 

where ip is a. conformal factor and jij a conformally related metric. We split the extrinsic 
curvature K^j into its trace K and a traceless part A^j according to 

Kij = Ai, + ^-ii,K = ip-^Aij + ^jijK. (11) 

In vacuum, the Hamiltonian constraint then becomes an equation for the conformal 
factor 

1 , ^ 1 1 
12 

Here Di and R are the covariant derivative and the Ricci scalar associated with the 
conformally related metric jij, and the covariant Laplace operator is = 'j^^DiDj. 

In the conformal thin-sandwich decomposition (see [321 [33] , as well as [211 [251 126] 
for reviews) we use the evolution equation for the spatial metric ([2]) to express A^^ as 

Here a = ip~^a and u^^ = dtjij, and the conformal Killing operator L is defined as 

[Lpy^ = D'(3^ + - '^fWkfi^. (14) 
The momentum constraint then becomes 

= {L(3y^ Dj In (a) + aDj (a-^u'^) + ^aip^D'K, (15) 



D^tp = -ipR + —ip^K^ - -i)-'^ AijA'^ . (12) 
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where = Dj{Lf3y^ is a vector Laplacian. In the extended conformal thin- 

sandwich formahsm we also combine the trace of the evolution equation for Kij with 
the Hamiltonian constraint to find an equation for the lapse a, 

(a^) = {^-ip-^A^.A^^ + ^^^K^ + ^i?) - ^''dtK + ^"^^'D^K. (16) 

The above equations form a set of equations for the lapse a, the shift (3^ and 
the conformal factor ■?/'. Before these equations can be solved, however, we have to 
make choices for the freely specifiable quantities for the conformal metric and its time 
derivative, 7^ and Uij = dt'jij, as well as for the mean curvature and its time derivative, 
K and dtK. 

To motivate our choices, consider a given stationary spacetime that, by virtue of 
being stationary, possesses a Killing vector. Let us choose an arbitrary spacehke slice 
through this spacetime. Projecting the spacetime metric gab onto this slice yields the 
spatial metric and, given the normal vector n", we can also compute the slice's 
extrinsic curvature Kij as well as its trace K. To assemble "conformal thin sandwich" 
data on this slice we also have to choose a time direction, which amounts to choosing a 
lapse and a shift. Aligning this time direction with the Killing vector, which yields the 
Killing lapse and shift, leads to a stationary slicing of the our stationary spacetime. We 
then have dt'jij = and dfK = 0. We also choose a conformal factor, so that we can 
compute the conformally related metric jij, and furthermore choose the conformal factor 
to be independent of time. This guarantees that Uij = 0. Clearly, then, our conformal 
factor i/j as well as the Killing lapse and shift must be solutions to the conformal thin- 
sandwich equations for the given ■jij and K, as well as Uij = and dtK = 0. Reversing 
the argument, we see that we can find this conformal factor and the Killing lapse and 
shift as solutions to the conformal thin-sandwich equations if we specify the given jij 
and K, as well as Uij = and dtK = as the freely specifiable quantities. 

To construct a stationary slice, or at least an approximately stationary slice through 
an approximately stationary spacetime, we should therefore choose Uij = and dtK = 
together with some and a K. Clearly, these choices cannot guarantee that the 
resulting slice will be stationary, since the conformal thin-sandwich formalism does not 
allow us to impose conditions on the first time derivative of the conformal factor, or the 
first time derivative of the tracefree part of the extrinsic curvature. However, the data 
is in some partial sense stationary because both dtjij = and dtK = 0. Choosing the 
mean curvature according to 

K=^, (17) 

in analogy to (iTl), further ensures dta = if we were to evolve these data using the 1 
+ log slicing condition We propose this condition as a "stationary l+log" slicing 
condition. In the following we will also assume conformal flatness, 'jij = rjij, where rjij is 
the flat metric in whatever coordinate system we choose. This conformal fiat assumption 
plays no real role in our result. 
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In addition to making choices for the freely specifiable quantities, we also need to 
impose suitable boundary conditions before the above equations can be solved. We will 
assume asymptotic flatness at infinity, which corresponds to 

lim V = 1 (18a) 
lim a =1. (18b) 

1 — >oo 

We will allow for a rotating coordinate system, so that the outer boundary condition on 
the shift becomes 

lim (3' = e'^^n.Xk. (18c) 

We construct black holes by excising the interior of coordinate spheres, and imposing 
the black hole equilibrium boundary conditions of [30] on the surfaces of these spheres. 
These conditions determine the conformal factor and shift on the black holes' horizons, 
but leave the lapse ojah undetermined. We discuss our choice for the horizon lapse 
below. 



3. Schwarzschild black holes 

Before exploring the slicing condition f|T7|) for rotating spacetimes in Section |H we 
will first examine non-rotating, spherically symmetric Schwarzschild spacetimes in this 
Section. 

To construct initial data, we solve equations (fT2l) . ( fTSl) and ( IT6l) under the 
assumption of spherical symmetry with Uij = = dtK, adopting the algorithm and 
numerical code described in [31]. The boundary conditions at spatial infinity are given 
by f ll8ap with Vtj = for nonrotating spacetimes. We also excise the black hole interior at 
a coordinate location tah- Following [30] we require the black hole to be in equilibrium, 
which results in boundary conditions on the conformal factor and the shift (3^; the 
lapse aAH, however, can be chosen arbitrarily. 

This situation is completely analogous to the maximal slices that we discussed at 
the end of Section [2^21 Had we imposed maximal slicing instead of the stationary 1+log 
slicing (fT7|) . then each value of oah would correspond to one particular parameter C in 
the family ([9]) of maximal slices of Schwarzschild. Since dynamical simulations settle 
down to the solution with C = 3v^M^/4, we would want to pick this particular solution 
by choosing a ah = 3v^/16. 

For the stationary slicing ( |T71) we see that we can again produce a family of solutions, 
parametrized by the horizon lapse oah- We again want to pick that member of this 
family to which the dynamical simulations settle down. This particular solution has 
been identified in [12], and we use this results to identify the corresponding value of 
ctAH- By inspection of their equation (9) we see that the lapse and the areal radius must 
simultaneously satisfy the equations 

3M/R -2 + 2a^ =0 (19a) 
2M/R - 1 + _ 2a = (196) 
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Figure 1. Graphs of the conformal factor ip^ the lapse a, and the shift for 
Schwarzschild in the stationary 1+log slicing (fT7|). 



at a certain critical radius for the solution to remain regular. Here M is the total mass 
of the spacetime, and, as before, R is the areal Schwarzschild radius. Solving these two 
equations yields ac = VTO — 3 and Rc = M/(4(vT^ — 3)). In spherical symmetry, the 
slicing condition ( fTTI) has a first integral 

= + (20) 

(see [15 )• Here C is a constant of integration, which plays the same role as the constant 
C in ([9]). In fact, this solution differs from ([9]) only in the exponential term. We now 
insert the above values of Oc and Rc into (1201) to find the value of C for the solution of 
interest 

C = J3-Vw)/2^2 ^ 1.24672M2. (21) 

16(v^- 3)3/2 ^ ^ 

Finally, we use this value together with R = 2M in fl20l) to obtain an equation for our 
value of the horizon lapse. Solving this equation iteratively yields 

aAH ~ 0.376179. (22) 

We adopt this value as a Dirichlet boundary condition for our differential equation. 
The solution of these equations is now completely determined, and we show numerical 
results for the conformal factor ip, the lapse a and the shift f3'^ in Fig. [H An alternative 
method for constructing these "trumpet" data is discussed in [H], who also provide an 
embedding diagram in their Fig. 2. 

In Appendix Appendix A we present some results from dynamical simulations of 
our Schwarzschild initial data and compare with evolutions of "puncture" data. 



The stationary 1+log slicing 



10 



4. Rotating spacetimes 

After exploring spherically symmetric spacetimes under the stationary 1+log condition 
f|T7|) . we next tried to construct solutions for binary black holes. Unfortunately, however, 
we were unable to do so. We have reason to believe that the problem is not of numerical 
nature - using maximal slicing K = instead of ( |T7I) . for example, our code reproduces 
the results of [30l [35] without problem - and instead believe that for rotating, non- 
axisymmetric spacetimes the condition (ITTl) does not allow any asymptotically flat 
solutions. We support this hypothesis with the following arguments. 



4.1. General considerations 

A standard approach to constructing 'quasi-stationary' initial data to the binary 
black hole problem goes as follows: One selects the independent data by specifying 
a background (conformal) metric and setting Uij = dt'jij = and dtK = 0. Because of 
the way the black holes orbit each other, one expects the best approximation to a Killing 
vector to be helical. Presumably, then, the time derivatives of jij and K are close to 
zero in a corotating coordinate system, which motivates our boundary conditions (118 ap . 
The orbital angular velocity Qj then appears in the equations as a free parameter that 
can be adjusted until a relativistic virial theorem is satisfied - namely the equality of 
the Komar mass with the ADM mass (see [36]) - indicating that the binary is in an 
approximately circular orbit. Most applications have adopted maximal slicing K = 0; 
here we consider the "stationary 1+log" slicing condition (fTTll instead. 



4-2. Spherically symmetry spacetimes 

Before proceeding, it is instructive to return to the spherically symmetric spacetimes of 
Section [31 

As above we will assume that the conformal factor falls off to one at infinity 
according to 

M , , 

^ = 1 + - + ---, (23) 

where M is the spacetime's ADM mass, and will analyze the asymptotic behavior of 
the lapse a and the shift (¥■ . Specifically, we will inspect the lapse equation (fT6l) . It 
is cleaner to assume that we have found a solution and use the conformal factor to 
transform to the solution frame. In this case the lapse equation becomes 

D'^a = aK'^Kij + (3'diK. (24) 

The solution of this equation can be written as a sum of a complementary solution acomp 
that satisfies the associated homogeneous equation, and a particular solution ctpart that 
"responds" to the source term. 

Ma Qa >^a /OK^ 

C^part 
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Here is a constant monopole moment, we assume that the dipole moment vanishes 
(which a suitable coordinate transformation will always accomplish), and we note that 
the quadrupole moment Qa has the standard flat-space harmonic function. We will 
assume that a falls off to one at infinity, i.e. we will assume n > 0. Inserting the ansatz 
( I25ll into (1241) we see that the left hand side scales according to 

- (26) 

where the symbol ~ indicates that asymptotically the leading-order term of the left 
hand side is proportional to the right hand side. 

With the assumption n > in (!25ll . the condition (|T7j) reduces to 

K ~ f3'dia, (27) 

Finally we assume that we can write the shift in the form 

= e'^%Xk + ^ + • ■ ■ (28) 

in Cartesian coordinates. Here the first term on the right-hand side is purely angular 
and accommodates the transformation from an inertial coordinate system into a rotating 
one. 

Under the assumption of spherical symmetry we may assume flj = in fl28l) . even 
though we will find below that we can always relax this assumption. We may also 
assume that the term in fl25]) is purely radial. In the following we abbreviate f3 = (3^ 
and B = B^. Inserting this into (p7|l yields 

^--.^ = ^.^ (29) 
where we have assumed > 1 for the moment, so that the leading term for a is given 
by the monopole term. To determine m we have to consider the momentum constraint 
( IT5|l . In spherical symmetry, this equation may be written as 

drO. 

r a ' ip 

(see equation (12b) in I34j)- Restricting attention to the leading order terms this reduces 
to 

^2/5 + ^ (^drf3 - ~ 0. (31) 

We can now insert the ansatz fl28|) and find m = 2 as an asymptotically flat solution. 
Inserting this into ( l29i) shows that K ~ r~^. We graphically verify that our numerical 
solutions of Section [3] satisfy this asymptotic fall-off behavior in Figure [2l From (fT3ll 
we also see that A^^ ~ (which is consistent with the observation that in spherical 
symmetry, traceless-transverse tensors are unique up to an overall scaling factor, fall off 
with 1/r^, and can be computed from {Lwy where is proportional to (l/r^,0,0)). 

Finally, we can determine n by inserting the decay rates of K, K'^^ and f3 into (12^ . 
The dominant term is the K'^^Kij term which falls off like 1/r^. This yields n = 4, which 
shows that our above assumption n > 1 was self-consistent for our purposes here. 





adrK (30) 
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Figure 2. The asymptotic behavior of the shift /3 and the mean curvature K in 
spherical symmetry. The sohd hnes denote the numerical results of Section [3l the 
dashed lines represent the asymptotic scalings (f35|) and (136]) . 



We could have obtained the scahng for a alternatively by observing that 

a'-(3' = l-—, (32) 

where, as before, R is the radial radius. If a approaches unity at infinity, and 
/3 ~ (9 we must have 

„=.l_2^ + 0f4) (33) 

or 



R Vi?^ 



-,A-^ + 0f4) (34) 



R \R^. 

If the spatial metric is in the standard Schwarzschild form, it is a straightforward 
calculation to show that D^{1 - 2M/Rf/'^ = 0. Therefore we have D'^a = O (r-^), 
exactly as required. 

In fact, we can also derive this asymptotic behavior analytically. Combining 
equations fl32l) and fl20l) we find 

CeV2 2.055M^ 

(35) 

We can now insert this into (fT7|) to find 

Ce^/^M^ 1.028M^ 

K ~ zrr ~ TTa • (36) 

2 R^ R'^ ^ ^ 

We include (1351) and (I36p in Figure [21 and find excellent agreement with our numerical 

results. 

Before closing this Section we point out that there was no need to assume flj = 
in (!28l) . The corresponding rotation term in the shift is purely angular, but since all 
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gradients must be purely radial in spherical symmetry, this term drops out both in fl27l) 
and ([23]). For example, if we orient the 2;-axis along Qi, then the rotational term in 
the shift is (3'^, but all derivatives along vanish in spherical symmetry. This means 
that we could, if we wanted, allow for a rotating coordinate system when constructing 
stationary 1+log slices of a spherically symmetric spacetime. This is a special case of a 
more general result for stationary axisymmetric systems which we will derive in Section 
14.41 For more general rotating systems, however, for which only a helical Killing vector 
exists, this does not seem to be possible. This is the content of the next subsection. 



4-3. Rotation 

We first observe that spacetimes describing rotating systems that are not axisymmetric 
generally do not possess a timelike Killing vector (e.g. they are radiating). Therefore we 
should generally not be able to find an exact stationary slicing of such spacetimes. 
Certain spacetimes of interest, e.g. binary black holes in approximately circular 
orbit, nevertheless possess an approximate helical Killing vector. For binaries, this 
helical Killing vector generates the orbit of the binary companions. If such an 
approximate helical Killing vector exists, we could apply the above formalism to 
construct approximately stationary slicings of such spacetimes. In Sections 12.21 and 
12.31 we identified the lapse and shift with the Killing lapse and shift, and based on 
these choices we set the time derivatives of the conformally related metric and the mean 
curvature to zero. This means that we assume the Killing lapse be of the form (l25ll 
and the Killing shift of the form fl28l) with Qj ^ 0. In practice, this angular velocity is 
determined as part of the iterative algorithm that is used to construct the solution. 

We again assume that we have found a solution, in which case the lapse a must 
satisfy equation (12^ . As before we assume that all the quantities belong to some 
weighted spaces, i.e., that each derivative introduces an extra power of r in falloff. It is 
easy to see that, at least generically, equations f|T7|) and flMj) are incompatible at large 
r. A simple counting argument goes as follows. Let us assume that the right hand side 
of (p^ falls off like l/r"^, say. Then (12^ implies that a must contain a term that decays 
like \jr^~'^ (the complementary part of a may contain terms that fall off even slower). 
Now insert this information in ( JT71) . The term dia decays like l/r"^"^, while /5, of course, 
diverges as r. Therefore K must decay like l/r"^"^. When this is inserted back into ([2 
we see that p^diK only decays like l/r"^~^, thus contradicting our original assumption 



X One might be concerned that the aKijK^^ and the (i^diK terms in ([24| could cancel each other, 
thereby avoiding the contradiction. We can show that this is impossible with the help of the following 
argument. Asymptotically, the shift /3* is dominated by /3' e^^^VljXk = Plou us therefore assume 
that aKijK^^ is canceled by Pl^^diK . But (31^^ is proportional to the i — 1, ra = Q magnetic vector 
spherical harmonic , which has axial parity. The angular components of diK have polar parity, 
and can be expanded into electric vector spherical harmonics . An integral of (il^^diK over any 
sphere must therefore vanish, indicating that Pl^^diK itself must take both positive and negative values 
(unless it is identically zero). The term aKijK^^ , however, is non-negative, meaning that it cannot be 
canceled by the (il^^diK term everywhere. 
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It is easy to see why the above argument does not hold for axisymmetric spacetimes. 
Let us assume, to leading order, that /5 is a rotation about the symmetry axis, which 
we align with the z-axis so that (3' ~ (0,0, /?<^) to leading order. But d^a = = d^K, 
since d^f, is a Killing vector. Therefore (3^dia and j3^diK vanish identically, rather than 
giving terms which decay slowly. In the presence of a true rotational Killing vector, 
therefore, equations f|T7|) and fl2^ are not incompatible, and it is possible to find 
asymptotically fiat, stationary l+log slices. This is not surprising, of course, since 
dynamical simulations of binary black hole mergers, using the 1+log slicing, settle down 
at late times to a time-independent and axisymmetric solution describing a single Kerr 
black hole, which clearly must satisfy both ( JTTll and ( |2^ . 

To account for the special case of axisymmetry we can revise the above argument 
as follows. Consider the nonaxisymmetric part of the right hand side of equation ( l24ll 
only. It is this part which we assume decays like This implies that a has a non- 

axisymmetric part which decays like l/r"^"^. This, when substituted into (fTTl) . implies 
that K has a nonaxisymmetric part which only decays like l/r"^"^. In turn, this means 
that the right hand side of equation (12^ can only decay like l/r"^"^, in contradiction 
to our above assumption. 



4-. 4- Stationary and axisymmetric spacetimes 

We now assume that the asymptotic behaviour of the spacetime is like that of the Kerr 
solution. In particular we assume that we have two linearly independent Killing vectors, 
namely a time translation Killing vector and a rotational Killing vector We 
will further assume that a slice E exists which is a stationary 1+log slice with respect 
to the timelike Killing vector T^, and that this slice is axisymmetric in the sense that 
the rotational Killing vector lies in the slice. We would then like to show that this slice 
is also a stationary 1+log slice with respect to any helical Killing vector 

H'^ = T^ + (37) 

where fl is now an arbitrary constant. 

We start by writing the Killing vectors in terms of their corresponding Killing lapses 
and shifts as 

= arn'' + (38) 

and 

^"k = a^n- + [31 (39) 

where n° is the normal on S. Since $^ is tangent to S, we must have = 0. Inserting 
( !38|) and (!39|) into (!37|) we can therefore write the helical Killing vector as 

Hl^ = aun'^ + Ph = c^m^ + Pt + ^Pl (40) 
and identify an = dr and = + 
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Given that S is a stationary 1+log slice with respect to we have 

K = (41) 
In order to show that S is also stationary 1+log with respect to H'^ we evaluate 

^ = = 2^r = + 2^^^^'^"^- ^^^^ 

Now consider a coordinate system with two unit vectors (ey)" and (e$)" that are aligned 
with the Killing vectors and In this coordinate system, the metric must then 
be independent of the coordinates T and $ associated with [exY and (e$)". Moreover, 
we have ut = {—g^^)~^^'^ ■ This implies that ax must be independent of $, and hence 
^'^daOtT = 0. The last term in equation ( H2i) therefore vanishes, so that the condition 
(H2|) holds if (HI!) does. This completes the proof (see also [37]). 

The above argument demonstrates both more formally and more generally an 
assertain that we made at the end of Section I4■2^ namely that there was no need to 
assume = in our analysis of spherically symmetric spacetimes. 



5. Summary and Discussion 

In this paper we propose and explore a "stationary 1+log" slicing condition, defined 
by equation ( fT7|) . for the construction of initial data. For stationary spacetimes, this 
condition automatically yields a stationary foliation when the data are evolved with the 
gauge conditions used in moving puncture simulations. The resulting slicing is time- 
independent and agrees with the slicing that is generated by dragging the initial data 
along the spacetime's time Killing vector. This is a desirable property, since it avoids 
pure "coordinate" evolutions that may otherwise introduce numerical error and noise, 
and that may contribute to the "spurious gravitational radiation" that is often observed 
at early times in binary black hole simulations. Given these considerations, this slicing 
condition would appear as an attractive alternative to maximal slicing, which is often 
used in the construction of binary black hole initial data. 

Unfortunately, it does not seem to be possible to construct initial data for binary 
black holes - or any binaries - in the stationary 1+log slicing. Spacetimes containing 
binaries in approximately circular orbit possess an approximate helical Killing vector 
that generates the binary's orbit. We provide an argument that shows that 1+log slices 
that are stationary with respect to such a helical Killing vector cannot be asymptotically 
flat, unless the spacetime possesses an additional axial Killing vector. As a consequence, 
the stationary 1+log slicing condition, without modification, seems useful only in the 
context of axisymmetry spacetimes, for example for rotating neutron stars. 

Our results also apply to generalizations of the 1+log slicing ([3]), for which the right 
hand side —2aK is replaced by —nf{a)K, where n is some arbitrary number and /(«) 
some non-zero and finite function of a (see |16j). 

It may be possible, however, to modify the stationary 1+log slicing condition in 
such a way that it can be applied to binaries as well. In particular, it seems promising 
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to apply 0171) only locally in a neighborhood of the binary, and match to maximal 
slicing K = asymptotically. Such an approach would still benefit from the slicing 
being stationary in the neighborhood of the binary, where the gravitational fields are 
strongest, but would avoid the problems associated with the asymptotic behavior of 
stationary 1+log slices. We plan to pursue this approach in the near future. 
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Appendix A. Dynamical evolution of Schwarzschild 

In this Appendix we present some results for dynamical evolution calculations for 
Schwarzschild spacetimes, starting both with "stationary 1+log" initial data ("trumpet 
data") and maximally sliced puncture data ("wormhole data"). While these results 
are not of direct relevance for the main arguments laid out in the paper, and while, 
being restricted to spherical symmetry, their scope is limited, they do demonstrate the 
potential of the stationary 1+log slicing condition to reduce numerical error. 

Dynamically evolving the initial data of Section[3]with the moving puncture method 
adopted in our code leads to one complication that is related to the black hole interior. 
Using the conformal thin-sandwich formalism for the construction of the initial data we 
have excised the interior by imposing boundary conditions on the black hole horizons. 
The initial data therefore cover only the black hole exterior. The moving puncture 
method, on the other hand, requires data everywhere, including the interior. One way 
to solve this problem is to fill the black hole interior with some arbitrary choices for 
the metric functions. As long as the constraint-violating data in the interior do not 
propagate to the exterior, this "junk" will not affect the outside of the black hole [38|[39]. 
at least in the continuum limit. In experiments for Schwarzschild black holes with the 
gauge conditions used in moving puncture evolutions, even the interior settles down to a 
constraint-satisfying solution that is determined by the exterior data, but not the initial 
choices inside the black hole. 

One source of error in this method arises from the junction of the exterior data to the 
arbitrary "junk" inside. In our finite-difference code, a grid stencil that is centered on a 
grid-point just outside the horizon may extend into the black hole interior. Numerically 
evaluated derivatives are therefore sensitive at these grid-points to the particular choices 
of the interior "junk" , and to how it is matched to the exterior solution (see [SHIEH]). The 
juncture of these data therefore results in a small numerical error, directly originating 
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Figure Al. Relative clianges in tlie ADM mass, AM = {M{t) - Mo)/Mo, wliere 
Mq is the initial ADM mass at t = 0, and where we compute the ADM mass by 
performing a surface integral on a sphere of coordinate radius r = 45.6Af. The solid 
line is the result of an evolution starting with puncture initial data representing a slice 
of constant Schwarzschild time ("wormhole" data), while the dashed line results from 
the evolution of our stationary 1+log initial data ("trumpet" data). 



from the black hole horizon, that can be controlled by improving the continuity of the 
data. In our implementation we construct the interior "junk" data using an 8th-order 
polynomial extrapolation of the exterior data into the interior. 

In Fig. lAll we compare the dynamical evolution of the normal puncture initial data, 
i.e. a slice of constant Schwarzschild time expressed in isotropic coordinates ("wormhole" 
data), with that of our data in stationary 1+log slicing ("trumpet" data). In both 
cases we use our three-dimensional code (i.e. we purposely do not take advantage of the 
spherical symmetry) that has been described in detail in jlOlHT]. As gauge conditions we 
use the 1+log slicing ([3]) together with the "unshifting shift" version of the "F^-freezing" 
condition 

dt(3' = (3/4)5^ (1.1a) 
dtB' = dtV ~ r]B' (1.16) 

with 7] = 1.0/M, where M is the black hole's ADM mass and F* = —dj'y^^. As initial 
conditions for these gauge conditions we choose 5* = as well as the Killing lapse 
and shift, as found from the solution of the conformal thin-sandwich equations. If is 
aligned with a Killing vector of the spacetime, then /3* = const is a solution to the shift 
condition (11. lap . 

Our code uses fourth-order spatial differencing and a fourth-order Runge-Kutta 
time stepping via the method of lines. We use the Carpet infrastructure ^] to 
implement the fixed mesh refinement (FMR) with nine levels of refinement. We resolve 
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the black holes at the finest level with a resolution of M/32, and our grid extends to 
256M at the outer boundary. 

As a measure of the numerical error we show in Fig. lAll relative changes in the 
ADM mass, evaluated as the surface integral 



on a sphere S of radius r = 45. 6M. Specifically, we graph AM = (M(t) — Mq)/Mq, 
where Mq is the initial ADM mass at t = 0. Strictly speaking, the ADM mass is 
defined only at spatial infinity. In numerical relativity codes whose grids do not extend 
to spatial infinity it is a common practice to evaluate an approximate value of the 
ADM mass at a finite separation instead. The resulting approximations converge to the 
correct ADM mass as the extraction radius approaches infinity; assuming reasonable 
fall-off properties this can be seen easily by using the Hamiltonian constraint to write 
the ADM as a volume integral. The errors shown in Fig. lAll are therefore, strictly 
speaking, not errors in the ADM mass, but close approximations and indicative of what 
would be observed in a dynamical simulation. 

Two observations are noteworthy. Firstly, adopting initial data on a slice of constant 
Schwarzschild time leads to deviations of the numerically determined ADM mass from 
its correct value that are larger than those resulting from the stationary l+log data. 
These errors are a consequence of the "wormhole" data collapsing into "trumpet" data. 
Secondly, we notice that these errors arrive at the extraction radius r = 45. 6M earlier 
than those for the stationary l+log data. The largest source of error for the "trumpet" 
data are those arising on the black hole horizon, as discussed above. As one might expect, 
these errors arrive at the extraction surface after a light-travel time of approximately 
t = 45M. For the "wormhole" data, on the other hand, the errors arrive earlier, because 
the transition from "wormhole" data to "trumpet" data occurs everywhere, including 
the exterior of the black hole (even though the size of these transitions decreases rapidly 
with the distance from the black hole). 

We also monitored the black hole's irreducible mass. While the relative deviations 
from the exact value were very small in both simulations (namely on the order of 10~^), 
the errors resulting from the "wormhole" initial data were about twice as large as those 
resulting from "trumpet" initial data. 

References 

[1] F. Pretorius. Evolution of binary black-hole spacetimes. Phys. Rev. Lett., 95:121101, 2005. 

[2] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and J. van Meter. Gravitational-wave extraction 

from an inspiraling configuration of merging black holes. Phys. Rev. Lett., 96:111102, 2006. 
[3] M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlochower. Accurate evolutions of orbiting 

black-hole binaries without excision. Phys. Rev. Lett., 96:111101, 2006. 
[4] M. Shibata and T. Nakamura. Evolution of three-dimensional gravitational waves: Harmonic 

slicing case. Phys. Rev. D, 52:5428, 1995. 
[5] T. W. Baumgarte and S. L. Shapiro. Numerical integration of Einstein's field equations. Phys. 




(1.2) 



Rev. D, 59:024007, 1999. 



The stationary 1+log slicing 



19 



[6] R. Beig and N. O Murchadha. Trapped Surfaces in Vacuum Spacetimes. Class. Quantum Grav., 
11:419, 1994. 

[7] R. Beig and N. O Murchadha. Vacuum spacetimes with future trapped surfaces. Class. Quantum 
Grav., 13:739, 1996. 

[8] S. Brandt and B. Briigmann. A simple construction of initial data for multiple black holes. Phys. 
Rev. Lett, 78:3606, 1997. 

[9] B. Briigmann. Binary black hole mergers in 3d numerical relativity. Int. J. Mod. Phys. D, 8:85, 
1999. 

[10] M. Alcubierre, W. Benger, B. Briigmann, G. Lanfermann, L. Ncrgcr, E. Scidcl, and R. Takahashi. 

3D Grazing Collision of Two Black Holes. Phys. Rev. Lett., 87:271103, 2001. 
[11] B. Briigmann, W. Tichy, and N. Jansen. Numerical Simulation of Orbiting Black Holes. Phys. 

Rev. Lett, 92:211101, 2004. 
[12] M. Hannam, S. Husa, D. PoUney, B. Bruegmann, and N. O'Murchadha. Geometry and Regularity 

of Moving Punctures. Phys. Rev. Lett., 99:241102, 2007. 
[13] M. Hannam, S. Husa, N. O Murchadha, B. Briigmann, J. A. Gonzalez, and U. Sperhake. Where 

do moving punctures go? 2006. 
[14] M. Hannam, S. Husa, F. Ohme, B. Briigmann, and N. O Murchadha. Wormholes and trumpets: 

the schwarzschild spacetime for the moving-puncture generation. 2008. 
[15] J. David Brown. Puncture evolution of schwarzschild black holes. 2007. 

[16] C. Bona, J. Masso, E. Seidel, and J. Stela. New Formalism for Numerical Relativity. Phys. Rev. 
Lett, 75:600-603, 1995. 

[17] M. Alcubierre, B. Briigmann, P. Diener, M. Koppitz, D. Pollney, E. Seidel, and R. Takahashi. 
Gauge conditions for long-term numerical black hole evolutions without excision. Phys. Rev. D, 
67:084023-+, 2003. 

[18] T. W. Baumgarte. Innermost stable circular orbit of binary black holes. Phys. Rev. D, 62:024018, 
2000. 

[19] K. A. Dennison, T. W. Baumgarte, and H. P. Pfeiffer. Approximate initial data for binary black 

holes. Phys. Rev. D, 74:064016, 2006. 
[20] H. P. Pfeiffer, D. A. Brown, L. E. Kidder, L. Lindblom, G. Lovelace, and M. A. Scheel. Reducing 

orbital eccentricity in binary black hole simulations. Class. Quantum, Grav., 24:59, 2007. 
[21] S. Husa, M. Hannam, J. A. Gonzalez, U. Sperhake, and B. Briigmann. Reducing eccentricity in 
black-hole binary evolutions with initial parameters from post-Newtonian inspiral. Phys. Rev. 
D, 77:044037, 2008. 

[22] R. Arnowitt, S. Deser, and C. W. Misner. The dynamics of general relativity. In L. Witten, 
editor. Gravitation: an Introduction to Current Research, page 227. Wiley, 1962. 

[23] R. M. Wald. General Relativity. The University of Chicago Press, Chicago, 1984. 

[24] G. B. Cook. Initial data for numerical relativity. Living Rev. Ret, 5:1, 2000. 

[25] T. W. Baumgarte and S. L. Shapiro. Numerical relativity and compact binaries. Phys. Rept, 
376:41, 2003. 

[26] E. Gourgoulhon. Construction of initial data for 3+1 numerical relativity. In M Alcubierre, editor. 
Proceedings of the VII Mexican School on Gravitation and Mathematical Physics. J. Phys.: Conf. 
Ser., in press. 

[27] F. Estabrook, H. Wahlquist, S. Christensen, B. Dewitt, L. Smarr, and E. Tsiang. Maximally 

sHcing a black hole. Phys. Rev. D, 7:2814 2817, 1973. 
[28] B. Reinhart. Maximal foliations of extended Schwarzschild space. J. Math. Phys., 14:719-719, 
1973. 

[29] R. Beig and N. O. Murchadha. Late time behavior of the maximal slicing of the Schwarzschild 

black hole. Phys. Rev. D, 57:4728-4737, 1998. 
[30] G. B. Cook and H. P. Pfeiffer. Excision boundary conditions for black-hole initial data. Phys. 

Rev. D, 70:104016, 2004. 

[31] T. W. Baumgarte and S. G. Naculich. Analj^tical representation of a black hole puncture solution. 



The stationary 1+log slicing 



20 



Phys. Rev. D, 75:067502, 2007. 
[32] J. W. York, Jr. Conformal 'thin-sandwich' data for the initial-value problem of general relativity. 

Phys. Rev. Lett., 82:1350, 1999. 
[33] H. P. Pfeiffer and J. W. York, Jr. Extrinsic curvature and the Einstein constraints. Phys. Rev. 

D, 67:044022, 2003. 

[34] K. Matera, T. W. Baumgarte, and E. Gourgoulhon. Shells around black holes: the effect of freely 

specifiable quantities in Einstein's constraint equations. Phys. Rev. D, 77:024049, 2008. 
[35] M. Caudill, G. B. Cook, J. D. Grigsby, and H. P. Pfeiffer. Circular orbits and spin in black-hole 

initial data. Phys. Rev. D, 74:064011, 2006. 
[36] E. Gourgoulhon and S. Bonazzola. A formulation of the virial theorem in general relativity. Class. 

Quantum Grav., 11:443, 1994. 
[37] N. O Murchadha. With commuting killing vectors, the lapse and shift of one killing vector are 

constants along the other. 2008. 
[38] Z. B. Etienne, J. A. Faber, Y. T. Liu, S. L. Shapiro, and T. W. Baumgarte. Filling the 

holes: Evolving excised binary black hole initial data with puncture techniques. Phys. Rev. 

D, 76:101503, 2007. 

[39] D. Brown, O. Sarbach, E. Schnetter, M. Tiglio, P. Diener, I. Hawke, and D. PoUney. Excision 

without excision. Phys. Rev. D, 76:081503, 2007. 
[40] J. A. Fabcr, T. W. Baumgarte, Z. B. Eticnnc, S. L. Shapiro, and K. Taniguchi. Rclativistic 

hydrodynamics in the presence of puncture black holes. Phys. Rev. D, 76:104021, 2007. 
[41] Z. B. Etienne, J. A. Faber, Y. T. Liu, S. L. Shapiro, K. Taniguchi, and T. W. Baumgarte. Fully 

general rclativistic simulations of black hole-neutron star mergers. Phys. Rev. D, 77:084002, 

2008. 

[42] www.carpetcode.org. 



